Magnetic orders and spin-flop transitions in the cobalt doped multiferroic 

Mni_^Co^W04 



Feng Ye,^'^ Songxue Chi,^ Jaime A. Fernandez-Baca,^'^ Huibo Cao,^ 
K.-C. Liang, ^ Yaqi Wang,^ Bernd Lorenz,^ and C. W. Chu^'^ 

^ Quantum Condensed Matter Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
'^Department of Physics and Astronomy, University of Kentucky, Lexington, Kentucky 40506, USA 
^Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996, USA 
^Department of Physics and TCSUH, University of Houston, Houston, Texas 77204, USA 
^Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 

(Dated: October 1, 2012) 

We present a comprehensive single-crystal neutron diffraction investigation of the Mni_xCoxW04 
with 0.02 < X < 0.30. At lower concentration x < 0.05, the system is quickly driven into the multi- 
ferroic phase with spin structure forming an elliptical spiral order similar to the parent compound. 
The reduction of electric polarization is ascribed to the tilting of the spiral plane. For x ^ 0.075, the 
magnetic structure undergoes a spin-flop transition that is characterized by a sudden rotation of the 
spin helix envelope into the ac plane. This spin structure persists for concentration up to x = 0.15, 
where additional competing magnetic orders appear at low temperature. For 0.17 < x < 0.30, the 
system experiences another spin- flop transition and recovers the low-x spiral spin configuration. A 
simple commensurate spin structure with ^= (0.5, 0, 0) is found to coexist with the incommensurate 
spiral order. The complex evolution of magnetic structure in Co doped MnW04 contrasts sharply 
with other transition metal ion doped Mni_xAxW04 (A=Zn, Mg, Fe) where the chemical substi- 
tutions stabilize only one type of magnetic structure. The rich phase diagram of Mni_xCoxW04 
results from the interplay between magnetic frustration and spin anisotropy of the Co ions. 

PACS numbers: 75.30.Kz,75.25.-j,61.05.F-,75.58.+t 



I. INTRODUCTION 



The observation of spontaneous electric polarization 
and magnetic control of ferroelectricity in perovskite 
manganite TbMnO^i has inspired much theoretical and 
experimental efforts searching for new magnetoelectric 
multiferroic materials due to their great technological 
and fundamental importance j^'i^ New magnetic multi- 
ferroics among transition metal oxides have since been 
discovered that include rare-earth manganite derivative 
i^MnOa and i?Mn205 {R: Y, Gd, Dy),^^ geometrically- 
frustrated triangular lattice CuTOs (T: Fe,Cr),^^io 
kagome lattice antiferromagnet Ni3V208.^^ In the con- 
ventional (type-I) ferroelectric materials, the electric po- 
lar state arises either from the covalent bonding between 
filled-shell oxygen atoms and the empty d-shell nonmag- 
netic transition metal ion {e.g., BaTiOs),^^ or the mver- 
sion symmetry breaking caused by the 65^ orbital (lone 
pair) that moves away from the centrosymmetric posi- 
tion {e.g., BiMnOs)^ Even if magnetic ions are present 
in these materials, the spins order at much lower temper- 
ature than the electric dipoles and the effect of magnetic 
transition on the dielectric constant is weak. In contrast, 
the new family of magnetoelectric multiferroics (type-II) 
exhibits an exceptionally strong sensitivity to an applied 
magnetic field that causes reversal and sudden flops of the 
electric polarization. ^^^^^ Such a level of control indicates 
the electric polarization is induced by the magnetic or- 
der, which typically has an incommensurate noncollinear 
spiral configuration. The onset of ferroelectricity corre- 
lates with the transition to the spiral spin order. The 



intimate link between these two order parameters marks 
the prominent and intriguing feature of the new class of 
multiferroics. 

The mineral hiibnerite MnW04 (monoclinic, P2/c) 
is one of the few multiferroics that is ideal for study- 
ing the complex spin orders caused by magnetic frus- 
trations and the interplay between magnetism and 
ferroelectricity^^^— Without chemical doping, the par- 
ent compound undergoes sequential magnetic transitions 
in zero magnetic field. ^^^^^ The system first enters a 
collinear spin structure (AF3) around 13.5 K with si- 
nusoidal modulation of the magnetic moment, the cor- 
responding incommensurate (ICM) wave vector appears 
at qs = (0.214,0.5,-0.457) and the moment of the Mn 
ions are confined in the ac plane at an angle of 35° to- 
wards the a axis. At 12.6 K, an ICM elliptical spiral 
spin order (AF2) sets in, accompanied by the sponta- 
neous electric polarization along the crystallographic b 
axis. The magnetic order at AF2 phase has the same 
magnetic wave vector q2 = qs while the helix spin struc- 
ture is characterized by the moment tilting out of the 
ac plane towards the b axis. As the system is further 
cooled below 7 K, the ICM magnetic order is replaced by 
a commensurate (CM) magnetic order (AFl) that also 
suppresses the electric polarization. The close proxim- 
ity of three different magnetic phases and the metamag- 
netic transitions prove the existence of significant mag- 
netic frustration in the system, as revealed by inelastic 
neutron scattering measurements. ^^^^^ Consistent with 
theoretical studies the investigation of spin wave 
excitations in the low-T ttii spin order indicates that 



the collinear spin state results from the balance of long 
range magnetic interactions between Mn^+ ions up to 
11*^ neighbors. Higher-order magnetic exchange inter- 
actions have strengths comparable to the nearest neigh- 
bor exchange coupling along the zigzag spin chain in the c 
axis and between the zigzag chains in the a axis direction. 
However, the magnetic interactions between the chains in 
the b axis are much weaker. This finding suggests that 
different magnetic structures are close in energy and com- 
pete for the magnetic ground state. Like other multifer- 
roic materials, the magnetic and ferroelectric phases in 
this system are expected to be tunable with small pertur- 
bations of magnetic field, ^'^'^^^ pressure, ^^^^ and even 
chemical substitution with various transition metal ions. 
MnW04 is known to form stable compounds when Mn is 
replaced by Fe,^^^^^ Co, Ni, Cu,^^ or Zn.^^ It was reported 
that the chemical doping with magnetic Fe ions stabilize 
the low-T collinear and CM spin order. ^^^^ Both neutron 
diffraction and electric polarization measurements show 
that the low-T magnetic ground state is completely con- 
verted into the collinear spin order (AFl phase) with 5% 
of Fe substitution, which is largely due to the increased 
Fe single ion anisotropy.^^ In contrast, substitution with 
nonmagnetic ions (Zn or Mg) has been shown to alter the 
magnetic ground state differently; replacement of only a 
few percent of nonmagnetic ions seems to be sufficient 
to suppress the commensurate AFl phase and stabilize 
the spiral spin order. ^^^^^ A similar effect was also ob- 
served in the Co-substituted compound Mni_xCoxW04 
by neutron diffraction and bulk magnetic measurements 
on polycrystalline samples. Song et ali^ found that the 
CM collinear state is replaced by the ICM spiral phase 
with 5% Co concentration; a spin-flop transition occurs 
with increasing Co doping with the corresponding spiral 
plane tilting away from the b axis. The modification of 
the magnetic structure is expected to switch the electric 
polarization direction from the b to the a axis as it was 
later confirmed in a single crystal of Mno.9Coo.1WO4.— »^ 
However, it was discovered that with further increasing 
of Co doping (~ 15%), the only measurable polarization 
was along the b axis.^^ This result indicates that a new 
type of spin structure responsible for the a axis polariza- 
tion is bound to a narrow range of Co substitution and 
the actual x — T phase diagram of the Mni_xCoxW04 is 
more complex than the one reported from powder diffrac- 
tion studies. 

To elucidate the nature of the magnetic ground 
state upon Co doping, we performed a comprehensive 
neutron single-crystal diffraction study of the doped 
Mni_xCoxW04. We observed a systematic evolution of 
the magnetic structure with increasing Co doping. The 
spin configurations for different phases are depicted in 
Fig. [H For lower concentration {x < 0.05), the system 
exhibits a spiral spin configuration similar to the undoped 
MnW04, but with decreased angle between the normal 
axis of the spiral plane and the c axis. At x ~ 0.075, 
the system undergoes a magnetic transition to a phase in 
which the spin helix flops into the ac plane. The ac spiral 
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FIG. 1. (Color online) (a) The crystal structure of the 
Mni_xCoxW04. The spin structures of (b) the collinear, 
commensurate (CM) AFl phase, spin moments are in the ac 
plane with an angle of ~ 35° to the a axis, (c) the collinear, 
incommensurate (ICM) AF3 phase, magnetic moments have 
the same direction as AFl but with modulated amplitude, 
(d) the noncollinear, ICM AF2 phase, one axis of the spiral 
ellipse lies in the ac plane, the other along the b axis, (e) the 
noncollinear ICM AF5 phase, the envelope of the spin helix 
lies in the ac plane, (f) the collinear, CM AF4 phase, spin mo- 
ments are in the ac plane with a angle of —50° to the a axis. 
Per rotation convention for the right-handed coordinates, the 
positive rotation angle with respect to the a axis is associated 
with the counterclockwise rotation when viewing the system 
along the negative direction of the axis of rotation. 



spin structure and the associated polarization Pa and Pc 
survive in the range of 0.075 < x < 0.15. For higher Co 
concentration, the system experiences a second spin-flop 
transition such that the spin order switches back to low- 
X spiral structure coexisting with a collinear CM mag- 
netic structure as observed in pure CoW04.^^ The com- 
plete phase diagram shown in Fig. [2] of Co-substituted 
MnW04 is obtained based on the neutron diffraction and 
bulk property measurements. The chemical substitution 
of the magnetic Co ion provides an unique method to fine 
tune the magnetic property that is capable of achieving 
magnetoelectric control with multiple value states. 

Single crystals of Mni_xCoxW04 with thirteen differ- 
ent compositions (0.02 < x < 0.3) have been grown in 
a floating zone optical furnace. We use powder x-ray 
diffraction to check the phase purity of the polycrys- 
talline feed rod before the crystal growth. No impurity 
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FIG. 2. (Color online) The phase diagram of Mni-xCoxW04 
as a function of temperature and Co concentration x. Due 
to the strong hysteresis of the magnetic order, transitions 
between different phases are identified with warming proto- 
col. Labels of various magnetic phases correspond to the spin 
structures displayed in Fig. 1. 
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FIG. 3. (Color online) Doping dependence of the lattice pa- 
rameter a,6,c and the angle /3 of Mni_xCoxW04 at room tem- 
perature characterized by x-ray powder diffraction. The crys- 
tal structure of all samples are refined with the monoclinic 
P2/c space group. Solid lines are guides to the eye. 



II. MAGNETIC ORDER AT LOW 
CONCENTRATION (0.02 <x < 0.05) 



phase could be detected within the resolution of the spec- 
tra. The chemical composition and the Co content of the 
single crystals were verified by energy-dispersive x-ray 
(EDX) measurements testing up to 15 different spots of 
a single crystal. As further confirmed by neutron diffrac- 
tion measurement, the refined Co content of all samples 
was close to the nominal composition. In the following 
sections, we will use the nominal composition to distin- 
guish between different substitution levels. Single crystal 
neutron diffraction experiments were carried out at the 
High Flux Isotope Reactor of the Oak Ridge National 
Laboratory. We used the HBIA, HBl, and HB3 triple 
axis spectrometers to study the doping and temperature 
evolution of the magnetic diffraction pattern. We chose 
an incident neutron beam with wavelength of 2.366 A 
and pyrolytic graphite (PG) crystals as monochroma- 
tor and analyzer. The crystals were aligned in several 
scattering planes to probe different magnetic reflections. 
For nuclear and spin structure determination, we used 
the HB3A four-circle diffractometer to collect both the 
nuclear and magnetic reflections with neutrons of wave- 
length 1.536 A at selected temperatures. The crystal 
and magnetic structure refinement were performed using 
the FULLPROF SUITE package. Magnetic representation 
analysis is performed to choose appropriate basis vec- 
tors to describe the various spin structures. The sample 
temperature was regulated either using a closed cycle re- 
frigerator (CCR) or liquid Helium cryostat. All samples 
under study belong to the monoclinic P2/c (No. 13) space 
group. With increasing Co concentration x, the lattice 
parameters of a, 6, c, and the angle of P systematically 
decrease over the range studied, as determined by the 
x-ray powder diffraction measurement shown in Fig. [3l 



We start with the neutron diffraction results of the 
lower Co concentration. The samples were aligned in 
a horizontal scattering plane defined by two orthogonal 
vectors of [1,0,-2] and [0,1,0]. Although the magnetic 
propagation wave vectors in the pure MnW04 is deter- 
mined to be (f = (0.214, 0.5, -0.457) that deviates slightly 
from the scattering plane used in the measurement, the 
coarse vertical resolution of the triple axis spectrometer 
is sufficient to capture the magnetic Bragg reflections and 
track their temperature and doping evolution. In the rest 
of the paper, we will use the wave vector (g, 0.5,— 2g), 
which is the projected value in the horizontal scattering 
plane, to denote the ICM magnetic Bragg peak unless 
specified otherwise. Figure |4] summarizes the thermal 
evolution of the ICM magnetic order at lower cobalt con- 
centration with X = 0.02,0.035,0.042, and 0.05. The in- 
tegrated intensities of magnetic scattering for these sam- 
ples grow monotonically below 13 K without any sign 
of anomaly at lower temperatures. Such behavior con- 
trasts with the undoped MnW04 where the multiferroic 
AF2 phase is sandwiched between the high-T AF3 and 
the low-T AFl phases, resulting in an abrupt suppression 
of the ICM magnetic scattering. Like the nonmagnetic 
Zn and Mg doping, ^^^^^ the multiferroic state associated 
with the noncollinear ICM phase is stabilized with Co 
substitution and becomes the ground state. The T de- 
pendence of the integrated intensities does not separate 
the phase boundary between the high-T AF3 and the 
multiferroic AF2 phases [see Fig. |4l^a)]. However, there 
are clear changes of the peak center across the phase 
transition that label the phase boundary between the 
collinear and noncollinear phases. In addition, the peak 
position is independent of the temperature in the AF2 
phase, but changes continuously in the AF3 phase, as 
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demonstrated in Figs. HJ^b) and HJ^c) for the x = 0.035 
and 0.05 samples. 
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FIG. 4. (Color online) (a) The T dependence of the ICM 
magnetic order parameters for Mni_xCoxW04 with x = 
0.02,0.035,0.042, and 0.05. The order parameter data are 
shifted vertically for clarity, (b) and (c) display the T de- 
pendence of the ICM peak centers upon cooling and warming 
for the X = 0.035 and x = 0.05 samples, respectively. Ar- 
rows indicate the transition between the high-T collinear and 
the low-T noncollinear phases. Contour plot of wave vector 
scans as a function of temperature for (d) x = 0.035 and (e) 
X = 0.05 upon cooling and warming. 



TABLE I. Refined parameters of the magnetic structure for 
the X = 0.02,0.035,0.042, and 0.05 samples in the non- 
collinear AF2 phase, mt and m±b denote the projected mo- 
ments along the two principle axes of the spin ellipse along 
and perpendicular to the b axis. is the angle between m±b 
and the a axis, or the angle between the normal vector n of 
the spiral plane and the c axis. The eccentricity of the spin 
helix is defined as e = (1 — rnl / rri\bY ^ to describe the devi- 
ation from a circular envelope. Notice the x — 0.075 sample 
undergoes rotation transition of the spiral plane on cooling 
and locks into the AF2 phase. Rp2 is the R factor calculated 
by Rp. = 100 E„(|G?b= - GLicI)/ E„ G?bs, where is the 
square of the structure factor for n observed reflections. 
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FIG. 5. Agreement of the magnetic structure refinements at 
(a) X = 0.02, (b) X = 0.035, (c) x = 0.042, and (d) x = 0.05. 
The calculated structure factors Fc^aic plotted against the 
experimentally observed Fobs- The corresponding spin struc- 
tures are shown in Fig. [TJd) and the inset of Fig. (6] 



The magnetic configurations of the AF2 phase for 
Mni_xCoxW04 at low x can be determined using sin- 
gle crystal magnetic diffraction experiments. As summa- 
rized in Table HI the refinement parameters indicate the 
spin structure is well described by the elliptical spiral as 
in the pure MnW04 and the doped Mni_xZnxW04 J^*^ 
Figure [5] shows the agreement between model calculation 
and experimental observations for all samples. Contrary 
to the Zn substitution, where the angle 9 between the 
normal vector n of the spiral plane and the crystallo- 
graphic c axis remains constant for Zn concentration up 
to X = 0.40, the value of shows a gradual decrease with 
increasing Co concentration [see Fig. El^a)]. The refine- 
ments also indicate a considerable deviation of the ellipse 
from a perfect circular envelope as quantified in the Ta- 
ble [H To correlate the spin structure and the ferroelec- 
tric properties, Fig. [6l^b) displays the T dependence of 
electric polarization P at different x. Clearly, a smooth 
decrease of the saturated polarization along the h axis is 
observed, a trend similar to the doping dependence of the 
angle 0. Within the microscopic picture of spin current 
or inverse Dzyaloshinskii-Moriya (DM) model, ^^^^ the 
expected electric polarization can be expressed as 



(1) 



where e^j is the unit vector connecting the neighboring 

spins Si and Sj^ and A is a constant related to the spin- 
orbit coupling and spin exchange interaction. ^^^^ The 
inset of Fig. [6] describes a general spin helix configuration; 
the normal vector n of the spiral plane has angle 6 with 
respect to the c axis and its projection in the ah plane has 
angle (j) towards the a axis. Using spherical coordinates, 
we can express n = sin cos (j)ex + sin sin (pey + cos Oez , 
where are the unit vectors along the i^^ Cartesian co- 
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ordinates {i = x^y^z). Since the angle between the crys- 
tallographic a and c axes is close to 90°, we will consider 
II a, e^y II 6, and || c for the sake of simplicity. The 
magnetic moments of two spins in the chemical unit cell 
for the helical state can be expressed as 

m{Ri^a) = my cos{27rqRi + + ^± sin(27r^^i + ^a)^ 

(2) 

where Ri^^ is the position vector of the Mn site a{= 1, 2) 
in the unit cell z, (f is the magnetic propagation wave vec- 
tor of the spiral structure, ^2 = ^ ^{Qz + 1), and my 
and m± are the long and short half axes of the magnetic 
ellipse, respectively. If we consider the dominant mag- 
netic exchange coupling is along the c axis zigzag chain, 
the electric polarization can be quantitatively written as 

Pi = Cimiimx sinyrg^ sin6>(— sin^e^^ + cos^e^y). (3) 

If, however, there is additional non- negligible exchange 
interactions between the chains along the a axis, the ex- 
tra contribution to the polarization is 
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P2 = C2Tn\\m j_ sin 27Tqx{— COS i 



sm U sm 



4>e.), (4) 



where Ci, C2 are constants independent of 6>, (j) and 
only related to the spin-orbit coupling and magnetic 
interaction, q^^ qz are components of the ICM mag- 
netic wave vector. The spin structure refinement in 
the AF2 phase reveals (j) = 180; this leads to a simpli- 
fied form of Pi = — Cim\\m± sin nqz sin Ocy and P2 = 
— C2Tn\\m± sin 27T qx cos Ocy^ both contributing to the b 
axis polarization. If there are no significant deviations 
of the ellipticity and the moment size (Table I), the mag- 
nitude of P depends solely on the angle 6. With both 
Pi and P2 taken into account, the total electric polar- 
ization becomes P5 ^ —m\\m±{Ci sin 6 + C2 cosO). It is 

evident that the polarization will change from Pi to P2 
upon variation of 6 from 90° to 0°. In addition, we note 
that IC2I < |Ci| due to a shorter exchange path along 
the zigzag chain direction (3.29 A along the c axis versus 
4.83 A along the a direction), the expression can be fur- 
ther simplified as P5 ~ —Cim\\m±sinO. One can then 
estimate that the saturated polarization Pij{x = 0.05) 
decreases to sin 15°/ sin 35° ~ 0A5Pb{x = 0) due to the 
rotation of the spiral plane, which is close to the mea- 
sured value. The close correlation between the rotation 
of the spin helix plane and the decreasing ferroelectric po- 
larization emphasizes that the doped MnW04 is indeed a 
prototypical multiferroic material with inverse DM mech- 
anism being the origin of the spontaneous electric polar- 
ization. 



III. MAGNETIC ORDER AT INTERMEDIATE 
CONCENTRATION (0.075 <x< 0.15) 

As the Co concentration increases to x = 0.075, fer- 
roelectric polarization measurement does not detect any 
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FIG. 6. (Color online) (a) The doping dependence of the 
angle 6 between spiral plane normal and the c axis for the 
Co doped MnW04. For comparison, the results of the Zn- 
doped samples are also shown. Inset shows the schematics 
of the elliptical spiral spin configuration in the AF2 phase. 
The normal vector n of the spiral plane has a angle of with 
respect to the c axis, the projection of ft on the ab plane 
has a angle of (j) to the a axis, (b) The T dependence of 
electric polarization Pb for x = 0,0.02,0.035,0.042 and 0.05 
Co-doping. 



significant b axis component within the resolution of the 
experiment. However, a large polarization was found 
along the a axis accompanied by a small component 
along c, which implies that the system undergoes ma- 
jor change in spin structure that causes the reorientation 
of the polarization. To correlate the polarization result, 
we studied the T dependence of the wave vector scans 
upon cooling and warming as displayed in Figs. [7l^a)-(b). 
One can clearly observe a strong hysteresis in the mag- 
netic scattering during thermal cycling. Upon cooling, 
the peak center shifts from qi ^ (0.234, 0.5, —0.468) that 
is associated with the high-T collinear AF3 phase and 
moves into a plateau with another ICM wave vector q2 ~ 
(0.232,0.5,-0.464) for 7 K< T <11 K [Figs. [T]; a) and 
[71(d)]. As the sample is further cooled to lower tempera- 
ture, the magnetic peak resumes shifting and finally locks 
into the low-T ICM phase with qs = (0.229, 0.5, -0.458). 
On warming, the same low-T magnetic order remains at 
the wave vector qs till T ~ 10 K and suddenly is con- 
verted to the high-T collinear AF3 phase. 

To understand the complex evolution of the magnetic 
order, the same sample was placed on the four-circle 
diffractometer to investigate the spin structure. We chose 
the cooling protocol and collected data at two character- 
istic temperatures T = 9 K and 4.5 K. Figures. [8ja) and 
mb) display the refined spin configurations. At T = 9 K, 
the magnetic spins form an ac spiral structure (termed 
as AF5 phase) with two principle axes of the ellipse ly- 
ing along the crystallographic a and c directions. The 
normal vector n of the new spiral plane is characterized 
hy = 90° and ^ = 90°. According to Eqn. ([1]), such 
a magnetic structure will cause the rotation of P5 into 
other directions. If we only consider that the dominant 
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FIG. 7. (Color online) The T dependence of the wave-vector 
scans of the ICM magnetic peak for x — 0.075 sample upon 
(a) cooling and (b) warming. The dashed lines in (a)-(b) 
are guides the eyes to illustrate the locked in peak position 
at the lowest temperature probed, (c)and (d) display the T 
dependence of the integrated intensity and the peak center of 
the ICM scattering. 



magnetic interactions are along the spin chain direction 
(c axis), the cross product of n and e^j would only induce 
the a axis electric polarization in a form of 



Pa = -Cirriamc^iniiqzex 



(5) 



where rria and rric are the projected moments of the spin 
helix along a and c. This prediction only agrees par- 
tially with the experimental observation where both Pa 
and Pc develop below 10 K.^^ To explain the presence of 
Pc7 we have to consider the interchain magnetic interac- 
tions along the a axis. Recent spin wave measurements 
in MnW04 have revealed that the interchain magnetic 
interactions are of the same order as the intrachain ex- 
change coupling in its collinear spin state. Although 
the results obtained are for the collinear spin structure, 
it is expected they will have similar strength when the 
system enters the noncollinear AF2 phase, based on what 
has been found in the pure CuFe02 and the multiferroic 
CuFei_xGax02i^i^ Taking the interchain magnetic in- 
teractions into account, the contribution to the electric 
polarization has the form of 



Pc = C2mamc sin 27: q^Cz 



(6) 



This prediction of small polarization along c is in ac- 
cordance with the bulk polarization measurement. We 
note that although the spins of the AF5 phase remain in 
the ac plane, the expressions of Eqns. [5] and [6] indicate 
a T-dependent eccentricity that can further modify the 
magnitude of the electric polarization. 
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FIG. 8. (Color online) The refined spin structure for x = 
0.075 sample at (a) 9 K and (b) 4.5 K. The corresponding 
agreement of the refinements are plotted on (c) and (d). Inset 
of panel (d) shows the agreement plot at T = 4.5 K using 
the same elliptical spiral at 9 K. The scattered data points 
indicate a poorer description of the spin configuration. 



With the sample cooled below 8 K, it was found that 
both Pa and Pc are suppressed. Such behavior suggests 
a gradual change of the spin structure that might be due 
to either the distortion of the elliptical envelope or the 
rotation of the spiral plane. Since both polarizations are 
proportional to maTric for an ac spiral structure, the in- 
creased eccentricity of the spiral ellipse could certainly 
yield a reduced bulk polarization. On the other hand, the 
rotation of the spiral plane away from the ac plane can 
also lead to the suppression of polarization along both 
directions. The spin structure determination with data 
collected at 4.5 K provides a definitive answer to sepa- 
rate those two possibilities. As shown in Fig. Efb), the 
system at base temperature shows a completely different 
spin structure from the one at 9 K. It possesses the con- 
figuration similar to the low Co concentration with the 
short axis of the spin ellipse along the b axis and the long 
axis residing in the ac plane. The canted angle between 
the spiral plane normal vector n and the c axis is ~ 6°. 
The result reveals that the x = 0.075 sample is located 
near the boundary between distinct spin configurations 
of the \ow-x spiral and the ac spiral structure. 

To understand the spin flop transition near Xc = 0.075, 
we recall that the spiral spin order allows the coupling 
linear in gradient of the magnetic order parameter (also 
known as Lifshitz invariant) with broken inversion sym- 
metry and induces a uniform electric polarization ac- 
cording to the Ginzburg-Landau approach. The energy 
gain originating from the nonlinear coupling between P 
and M (- 7P • [M(V • M) - (M • V)M]) and the corre- 
sponding electric polarization is proportional to the angle 
0. With increasing Co concentration, both the electric 
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FIG. 9. (Color online) The T dependence of the wave- vector 
scans of the ICM magnetic peak for x = 0.10 sample upon 
(a) cooling and (b) warming. Notice the weak commensurate 
AFl reflections at g = (0.25,0.5,-0.5) is present near 11 K 
for both cooling and warming protocols, (c) and (d) display 
the T dependence of the integrated intensity and the peak 
center of the ICM scattering. 



polarization P and the energy gain decrease because of 
the decreasing angle. It is not surprising that exists a 
critical concentration Xc, around which the system can 
no longer gain energy to maintain the magnetic struc- 
ture. Such instability will then bring about a spin flop 
transition to the observed ac spiral which can reduce the 
free energy further because of the sizable polarization. 
Compared to other multiferroics, the continuous change 
in the magnetic wave vector from the transition between 
different spin configurations of the Co doped MnW04 
appears rather rare. The gradual transition highlights 
that the high-T spiral and the low-T ac spiral spin struc- 
tures are nearly degenerate in energy such that a small 
change in temperature would drive the system between 
competing magnetic states characterized by a continuous 
rotation of the spiral plane. 

Figures. [9l^a)-|9l^d) show the neutron diffraction results 
on the X = 0.10 sample. The hysteresis of the magnetic 
order upon thermal cycling is still discernible. However, 
the shift of the magnetic wave vector is much smoother 
upon cooling and does not show the plateau as seen at 
X = 0.075. Accordingly, electric polarizations show con- 
tinuous growth upon cooling. The saturated Pa and Pc 
reach 100 and 30 /iC/m^ at 4 K, respectively. Pa exceed 
the maximum 6-axis polarization of MnW04 by nearly 
a factor of two. This result is consistent with the pre- 
dicted P in Eqn. (0 assuming unchanged coefficient Ci 
in the AF2 and AF5 phases. On warming, the low-T spin 
structure remains till 11 K and a transition to the high- 
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FIG. 10. (Color online) The T dependence of the wave- vector 
scans of the magnetic peaks for x = 0.12 sample upon (a) 
cooling and (b) warming. The system undergoes ICM AF3, 
CM AFl, and ICM AF5 phases as cooling. Similar scans on 
the X = 0.135 sample are presented on (c) and (d). Only CM 
AFl and ICM AF5 orders appear as the major phases. In 
(c), a minor AF2 phase with q = (0.23, 0.5, —0.46) appears at 
5.8 K < T < 7.6 K. 



T collinear AF3 phase occurs. The lack of an abrupt 
change in the ICM wave vector suggests the modification 
of the spin structure is not as drastic as the x = 0.075 
sample and the same ICM spin structure is preserved to 
the lowest temperature. In a recent neutron diffraction 
study on the x = 0.10 sample, Urcelay-Olabarria et al. 
obtained a similar result and found that the magnetic 
spins remain in the ac plane at all temperatures while 
the eccentricity reduces from e = 0.66 at 9 K to 0.42 at 
2 K, suggesting the spiral ellipse becomes more circular 
at lower temperature.^^ This decrease of the eccentricity 
upon cooling further increases the electric polarization, 
as expressed in Eqns. (|5]) and ([6]). Another noteworthy 
feature in the x = 0.10 sample is that the energy of the 
ac spiral structure is close to other competing state in- 
cluding the collinear AFl phase. This can be appreciated 
by the weak magnetic reflection at the CM wave vector 
q = (0.25,0.5,-0.5) near 11 K [see Figs. [H a) andlljb)]. 

The evolution to other competing magnetic orders is 
manifested at higher Co concentration. Figure [10] dis- 
plays the T dependence of the wave- vector scans for two 
samples with x = 0.12 and x = 0.135. At x = 0.12, the 
first magnetic order present in the scattering plane con- 
sisting of the [1,0,-2] and the [0,1,0] directions is the 
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FIG. 11. (Color online) The T dependence of the wave- vector 
scan near the magnetic orders of the x — 0.15 sample along 
(a) the [if,0.5, -2if] and (b) the [if, 0, 0] directions, (c)-(f) 
are the T dependence of the integrated intensities for the AF2, 
API, AF4, and AF5 phases. 



collinear AF3 phase, which locks into the commensurate 
AFl phase at lower temperature. With the sample cooled 
below 9.0 K, the previously observed AF5 phase sets in 
and extends to the lowest temperature without any no- 
ticeable shift of the magnetic wave vector. There is a nar- 
row temperature window in which both the AFl and AF5 
phases are present between 6.8 K < T < 9.0 K. On warm- 
ing, the reverse order of the magnetic phases is observed 
except the coexisting region shifts up to 7.5 K < T < 
9.3 K. Interestingly, the order of appearance for the com- 
mensurate AFl and the multiferroic AF5 phases is ex- 
actly opposite to the pure MnW04, where the AFl phase 
occurs at lower temperature. At x = 0.135, the transition 
between the collinear AF3 and the AFl phases vanishes. 
There is very weak scattering near the ICM wave vec- 
tor q = (0.225,0.5,-0.45) that exists in a very narrow 
range of 10.7 K< T <11.3 K, which was later identified 
as the AF2 phase. The strongest intensity is only one 
percent of that for the AFl phase and is too weak to 
induce any detectable polarization signal. We also ob- 
served the coexistence of the multiferroic AF5 phase and 
collinear AFl phase. The coexisting region downshifts 
to 5.5 K< T < 7.2 K for cooling and 7.0 K< T <7.8 K 
for warming. At lower temperature, the AFl phase is 
completely suppressed. 

The X = 0.15 sample is probably the most complex sys- 
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FIG. 12. (Color online) The reciprocal space mapping near 
the CM and ICM magnetic reflections at (a) 2, (b) 4, (c) 6.5, 
and (d) 10 K. At least four distinct magnetic reflections are 
observed in this scattering plane at low temperatures. 



tem with a minimum of five coexisting magnetic phases 4^ 
Figure [11] summarizes the thermal evolution of the vari- 
ous magnetic orders probed at two scattering planes. Fig- 
ure fTU a) shows the T dependence of wave- vector scans 
along the [i^, 0.5, -2i^] direction where the AFl, AF2, 
and AF5 can be surveyed. Figure [TTJb) displays the 
scans along the [H, 0, 0] direction where the commen- 
surate AF4 magnetic order with q = (0.5,0,0) can be 
examined. Upon cooling, the AF4 phase first appears 
around 17 K and the intensity increases continuously till 
10 K. A sharp drop in its intensity is accompanied by the 
simultaneous development of the commensurate AFl and 
multiferroic AF2 phases. With the sample cooled below 
6.6 K, the CM AFl phase is also suppressed and an ad- 
ditional ICM AF5 phase develops at lower temperature. 
Figures. [mc)-[TT]ff) summarize the T dependence of inte- 
grated intensities for the four major magnetic (the com- 
mensurate AFl, AF4 and incommensurate AF2, AF5) 
phases in the x = 0.15 sample. Since the magnetic in- 
tensity of the AF2 phase closely follows the b axis po- 
larization P5, it is speculated that this state has a spiral 
spin structure similar to the x < 0.05 samples. Notice 
even at the lowest temperature, there are finite magnetic 
scattering from the commensurate AFl and AF4 phases. 

To confirm the coexistence of various magnetic orders, 
we have performed an extensive survey in the recipro- 
cal space at selected temperatures, 2.0, 4.0, 6.5, and 
10 K, as shown in Fig. [121 The appearance and dis- 
appearance of competing magnetic orders is evident. At 
10 K, only the AF2 with q2 = (0.225,0.5,-0.45) and 
AFl with qi = (0.25,0.5, —0.5) phases are present, while 
more magnetic Bragg peaks appear at lower temperature. 
A distinct magnetic reflection appearing as the shoul- 
der of the AF5 Bragg peak becomes visible in the low-T 
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FIG. 13. (Color online) The T dependence of the integrated 
intensity for (a) the AFl phase with q= (0.25, 0.5, —0.5) and 
(b) the AF4 phase with q= (0.5, 0, 0) at x = 0.12, 0.135, and 
0.15. Both the AFl and AF4 states survive at low tempera- 
ture for the X — 0.15 sample. 



mapping [see Figs. [I2fa) and [l2fb)]. The peak is lo- 
cated between the AF2 and AF5 phases with wave vector 
q = (0.23, 0.5, —0.46). The exact nature of this magnetic 
order remains unknown since its reflection is too close to 
the neighboring AF2 and AF5 Bragg peaks. This un- 
known magnetic phase is already present at x = 0.135, 
but only exists within a narrow temperature range [see 
Fig. [TOT c)]. This magnetic fluctuation continues to grow 
in intensity upon cooling at x = 0.15. Together with the 
already identified commensurate AF1,AF4 phases and in- 
commensurate AF2,AF5 phases, there are five magnetic 
orders at the lowest temperature. Such remarkable co- 
existence of many competing magnetic phases marks the 
X = 0.15 sample as the most frustrated system. 

Table II lists the refined magnetic structures for the 
X = 0.075, 0.10, 0.12, and 0.135 samples in the AF5 phase 
as well as the commensurate AFl and AF4 phases for 
the X = 0.135 sample. As mentioned in Sec. IIIH the 
X = 0.075 sample is located near the phase boundary be- 
tween the \ow-x spiral and the ac spiral structures; its 
magnetic structure can be refined as an ac spiral con- 
figuration only for T > 8 K. In contrast, the magnetic 
ground states of the x = 0.10,0.12 samples are well de- 
scribed by the same ac spiral structure at low temper- 
atures. Although the ICM magnetic order seems to be 
the only low-T phase for the x = 0.135 sample, it cannot 
be refined by a pure ac spiral structure implying the de- 
viation from that configuration. To get good agreement 
between the observation and model calculation, a com- 
bined low-x spiral and ac spiral structure is chosen to fit 
the experimental data and yields satisfactory result. The 
spin structure is best characterized as a modified ac he- 
lical structure that the normal vector of the spiral plane 
tilts way from the b axis. Such spin order results in a 
reduction of electric polarization and is consistent with 
the bulk measurement that both Pa and Pc decrease for 
X > 0.10.ii 



TABLE II. Refined parameters of the magnetic structures at 
X = 0.075,0.10,0.12, and 0.135. The spin configurations of 
the X = 0.075 sample at 9 K and the x = 0.10, 0.12 samples at 
5 K can be refined as an ac spiral structure with its principle 
axes aligning along the a and c directions. The x = 0.135 
sample is refined to be the AF4 phase at 12 K, the AFl phase 
at 9 K and the modified AF5 phase at 5 K. The moment 
direction has a angle of —33° with respect to the a axis in the 
AFl phase and —52° to the a axis in the AF4 phase. At the 
AF5 phase, the spin helix has one of its two principle axes in 
the ac plane and the other in the be plane. 
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With Co concentration x > 0.10, we have observed the 
expansion of both commensurate AFl and AF4 phases 
as shown in Fig. [131 Overall, samples in this doping 
region form the commensurate AF4 spin structure with 
q4 = (0.5,0,0) at higher temperature and enter directly 
the collinear AFl phase upon cooling, which is different 
from pure MnW04. With increasing x, the transition to 
the AFl state moves to lower temperature while the tran- 
sition to the AF4 phase shifts to higher one. Although 
both phases at x = 0.10,0.12,0.135 are completely sup- 
pressed at low-T, they survive for the x = 0.15 sample 
indicating the collinear spin structures gradually become 
the stable magnetic ground state at large x. Magnetic 
structure refinements in this doping range reveal that the 
spin moments in the AF4 phase are confined in the ac 
plane, with an angle of —50° towards the a axis (see Ta- 
ble II). This spin reorientation is again due to the strong 
anisotropy of Co^+ ions^^^^^ that locks the Mn^+ spins 
in the same direction and makes the collinear spin struc- 
ture more favorable with increasing x. The spins in the 
collinear AFl phase have a angle of —33° with respect to 
the a axis, which is different from pure MnW04, and is 
probably due to the pinning of the high-T AF4 magnetic 
structure. 



IV. MAGNETIC ORDER AT HIGH 
CONCENTRATION (0.17 <x< 0.30) 

Finally, we focus on the magnetic structures for x > 
0.17. The bulk polarization measurements show no de- 
tectable Pa and Pc. Instead, the polarization is pointing 
to the b axis, the same direction as observed at lower 
Co concentration. Thus, the switch of the polarization 
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FIG. 14. (Color online) (a) The T dependence of the wave- 
vector scans of the magnetic orders along the {H,0, —2H) 
direction for x = 0.17, 0.20, 0.25, and 0.30. 



suggests another major modification of the spin struc- 
ture. Figure [T4l compares the T dependence of the wave- 
vector scans across the ICM peak of four Co doped sam- 
ples with X = 0.17,0.20,0.25, and 0.30. The scattering 
profile in this doping regime exhibits different charac- 
ter. Unlike the coexistence of various competing mag- 
netic orders in the intermediate doping regime, there is 
only one ICM magnetic reflection at the wave vector of 
q ^ (0.22,0.5,-0.44). The transition temperature de- 
creases from 9.6 K at x = 0.17 to 4.4 K at x = 0.30. The 
scattering intensity of the ICM magnetic order is also 
suppressed with increasing x as displayed in Fig. [TST a). 
Further survey in reciprocal space reveals one strong 
collinear AF4 phase that is established at higher tem- 
peratures and persists to the lowest temperatures [see 
Fig. [TST b)]. The transition temperature increases with 
Co concentration and reaches 25 K for the x = 0.30 sam- 
ple. For all samples studied, the magnetic intensities 
of the AF4 phase exhibit a kink at temperatures corre- 
sponding to the onset of the mentioned ICM magnetic 
order indicating the competition between the CM and 
ICM phases. 

Since all samples in this doping region show similar 
magnetic properties except the transition temperature, 
we chose the x = 0.17 sample for the crystal and mag- 
netic structure refinement and expect the other doped 
samples have smooth evolution of the spin structure. As 
demonstrated in Figs. [14] and [151 there are two major 
magnetic phases with ICM and CM wave vectors for 
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FIG. 15. (Color online) (a) The T dependence of the inte- 
grated intensities of the AF2 phase for x = 0.17,0.20,0.25, 
and 0.30. (b) The T dependence of the integrated intensities 
of the AF4 phase with q = (0.5, 0, 0) for the same samples. 
Arrows label the transitions where the low-T AF2 phases set 
in. The magnetic scattering intensities are normalized to the 
intensities from the nuclear reflections for comparison. 
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FIG. 16. (Color online) The spin configurations of (a) the low- 
T incommensurate AF2 phase, (b) the high-T incommensu- 
rate AF4 phase, and (c) the conical spin order resulting from 
the superposition of the AF2 and AF4 phases. Agreement 
plots for (d) the incommensurate AF2 phase at 5 K, (e) the 
commensurate AF4 phase at 5 K, and (f) the commensurate 
AF4 phase at 11 K. 



X = 0.17. We collected 193 nuclear reflections at 5 K 
for the structural determination. One set of magnetic 
reflections for the low-T ICM magnetic structure were 
collected to refine the spin structure and two sets of 
magnetic reflections with q = (0.5,0,0) were collected 
at 5 and 11 K separately to investigate how the collinear 
AF4 phase is affected by the low-T ICM magnetic or- 
der. Figures. [TOT a) and [TOT b) show the corresponding 
magnetic structures of the ICM spiral order and the CM 
AF4 phase. At 11 K, only the collinear AF4 phase ex- 



11 



TABLE III. Magnetic structure parameters at x = 0.17. The 
sample has colhnear AF4 phase with q = (0.5,0,0) at 11 K 
and 5 K. At T = 5 K, additional AF2 spiral structure is 
formed. rUa and rric are the spin components along the a and 
c axes in the collinear phase. The real m±b and imaginary irib 
denote the projected moments of the long and short spiral 
ellipse lying in the ac plane and along the b axis. is the 
angle of the m±b with respect to the a axis. 
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ists, the magnetic spins form in a configuration identical 
to C0WO4, where the moments lie in the ac-plane, with 
an angle of ^ ~ —50° towards the a axis. With the sam- 
ple cooled below 10 K where the ICM order sets in, the 
magnetic structure of the collinear spin order is not mod- 
ified; the spins remain in the same direction and the total 
moment of 2.42(4) /i^/site at 5 K is almost the same as 
2.39(3) /i^/site at 11 K (see Table III). On the other 
hand, attempting to refine the low-T ICM order using 
the ac spiral structure gives a poor fit to the collected 
data. Consequently, we adopted the spin configuration 
at low Co concentration because such spin structure pro- 
vides an electric polarization along the b axis. As shown 
in Fig. [TCT d), this model provides a good description of 
the experimental data. The spin moment of the spiral 
state varies from 2.36 /iB to 2.92 /i^, which is comparable 
with the moment size at the AF4 phase. The observation 
of spiral order with the helix plane similar to the \ow-x 
case implies that the x = 0.17 sample is located at an- 
other phase boundary where the system undergoes a sec- 
ond spin flop transition consistent with P || However, 
there is one apparent difference between the x < 0.05 
and X > 0.17 samples. There is only one ICM AF2 phase 
in the low Co concentration samples, while the ICM AF2 
phase appearing in the high-x samples coexists with the 
commensurate AF4 phase that is established at higher 
temperature. The simultaneous presence of two magnetic 
phases is similar to the colossal magnetoresistance related 
manganite PrcyCacsMnOa (PCMO), in which both the 
ferromagnetic and antiferromagnetic components are ob- 
served at low temperature^ This can be interpreted ei- 
ther as a canted antiferromagnetic structure or coexisting 
FM and AFM phases. Similarly, neutron diffraction data 
alone can not differentiate whether the observed coexis- 
tence of AF2 and AF4 orders in Mno.83Co.i7W04 arises 
from two separated phases, each with distinct magnetic 
wave vector; or if they originate from one single phase 
with two- A: magnetic structure. Definitive identification 
would require a spatially sensitive probe in conjunction 
with the electric polarization measurement. Our neu- 
tron diffraction result at x = 0.17 is consistent with 
a recent study on the x = 0.20 sample, where the co- 



existing collinear AF4 and multiferroic AF2 phases are 
revealed. Those authors concluded that the superposi- 
tion of the competing AF4 and AF2 magnetic structures 
leads to a conical antiferromagnetic order that is depicted 
in Fig. [TBTc). 
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FIG. 17. (Color online) The doping dependence of the incom- 
mensurability of the low-T magnetic structures. The low-T 
spin orders are labeled with different symbols. The shadow 
region for the x = 0.075 and x = 0.10 samples indicates the 
samples undergo continuous change of the magnetic wave vec- 
tors upon cooling. Solid lines are guides to the eye. For 
X — 0.15 and 0.17, there are additional magnetic Bragg re- 
flection between the known AF2 and AF5 phases. 

The comprehensive single-crystal neutron diffraction 
measurement, in combination with the magnetic prop- 
erty and polarization measurements, make it possible to 
construct the phase diagram of the Mni_xCoxW04 as a 
function of Co concentration and temperature. Unlike 
other transition metal ion doped MnW04, where only 
one type of spin configuration is stabilized, a rich va- 
riety of spin structures and complex evolution between 
different phases are observed as the Co concentration is 
increased. The spin anisotropy of the Co ions plays a 
vital role in defining the low-T magnetic structures. Al- 
though confined in the same ac plane, the spin easy axis 
in C0WO4 is —45° from the a axis, and is nearly 90° 
away from the easy axis direction in MnW04.^^'^^ With 
increasing x, the long axis of the spiral ellipse that ini- 
tially has a positive angle towards the a axis will tilt 
gradually to the negative direction due to the single ion 
anisotropy of the Co^+. The rotation of the spin helix 
plane leads to a decrease of the electric polarization |P| 
that is compatible with the magnetic symmetry, as well 
as the energy gain that is coupled to the |P|. At the 
critical concentration oi x = 0.075, the system cannot 
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gain enough energy to maintain the multiferroic phase, 
thus results in the spin flop transition. The new ac spiral 
structure helps the system lower the free energy because 
of the large value of the ferroelectric polarization. For 
0.075 < X < 0.15, the presence of the Pa and Pc instead 
of Pb is consistent with an ac spiral spin structure, and 
highlights the significant intrachain as well as the inter- 
chain interactions. Although the samples in the interme- 
diate doping exhibit similar spiral structure, the maxima 
of the electric polarizations only occur near x = 0.10. 
It results mainly from the rotation of the spiral plane 
as exemplified by the neutron diffraction data from the 
X = 0.075,0.135 samples. The perfect ac spiral order is 
realized at x ^ 0.10 that induces the largest polariza- 
tion. Further increasing x not only causes the deviation 
from a pure ac spiral configuration, but also introduces 
collinear AFl and AF4 spin orders that reduce the ef- 
fective moment of the spiral structure responsible for the 
polarization. A second spin-flop transition takes place 
with X > 0.17 and leads to a similar spiral structure as 
in the low Co concentration. This phase coexists with a 
CM collinear AF4 phase with spin configuration similar 
to C0WO4. The gradual suppression of the electric po- 
larization in this doping region is mainly caused by the 
increasing collinear AF4 order, other than the rotation 
of the spiral plane, as observed at x < 0.05. 

The phase diagram is characterized by three well- 
defined regions distinguished by different spin spirals. 
However, the evolution of the magnetic spin structure 
within an individual region is gradual. This can be better 
appreciated by examining the concentration dependence 
of the incommensurability of the low-T noncollinear or- 
der, as shown in Fig.[T3 The magnetic wave vector of the 
ICM order does not exhibit a lock-in value over the wide 
range of Co concentration. Instead, it varies smoothly 



within each of the three regions. For the samples near 
the X = 0.075 phase boundary, the ICM magnetic struc- 
ture is unstable such that a small temperature variation 
will cause the rotation as well as the periodicity change of 
the spiral structure. All these observations reinforce that 
the magnetic structure results from the delicate balance 
between the competing exchange interactions and spin 
anisotropy of the transition metal ions. 

Without chemical substitution, MnW04 appears to be 
a highly frustrated system where the magnetic and fer- 
roelectric properties can be modified by external stimuli 
like the magnetic and electric field. ^^^^ The introduction 
of Co ions with distinct spin anisotropy provides another 
way to fine tune the magnetic ground states. This is sim- 
ilar to the rare-earth multiferroic manganite i^MnOs, in 
which the ferroelectric polarization is enhanced by the 
magnetic order of rare earth elements, except that the 
tuning parameter is on the same magnetic site in the case 
of Mni_xCoxW04. We hope the current experimental 
study will inspire further theoretical effort to understand 
the magnetic and ferroelectric order parameters in this 
doped system. Most importantly, such work will provide 
a new pathway to design and synthesize magnetoelectric- 
control materials with multiple magnetic and ferroelec- 
tric ground states. 
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